Early onset of ground-state deformation in the neutron-deficient polonium isotopes 



T.E. Cocolios 1 - 2 , W. Dexters 1 , M.D. Seliverstov 1 ' 3 - 4 , A.N. Andreyev 1 - 5 , S. Antalic 6 , A.E. Barzakh 3 , 

B. Bastin 1 *, J. Buscher 1 , I.G. Darby 1 , D.V. Fedorov 3 , V.N. Fedosseyev 7 , K.T. Flanagan 8 ' 9 , 

S. Franchoo 10 , S. Fritzsche 11 ' 12 , G. Huber 4 , M. Huyse 1 , M. Keupers 1 , U. Koster 13 , Yu. Kudryavtsev 1 , 

E. Mane 87 , B.A. Marsh 7 , P.L. Molkanov 3 , R.D. Page 14 , A.M. Sjoedin 7 - 15 , I. Stefan 10 , J. Van de 

Walle 1 ' 2 *, P. Van Duppen 1 , M. Venhart 1 - 16 , S.G. Zemlyanoy 17 , and M. Bender 18 , P.-H. Heenen 19 
Instituut voor Kern- en Stralmgsfysica, K. U. Leuven, B-3001 Leuven, Belgium 
2 ISOLDE, CERN, OH- 1211 Geneva 23, Switzerland 
3 Petersburg Nuclear Physics Institute, 188350 Gatchma, Russia 
4 Institut fur Physik, Johannes Gutenberg Universitat, D-55099 Mainz, Germany 
5 School of Engineering and Science, University of West Scotland, Paisley, 
PA1 2BE, United Kingdom and Scottish Universities Physics Alliance (SUPA) 
6 Department of Physics and Biophysics, Comenius University, Bratislava 842 48, Slovakia 

7 EN Department, CERN, CH-1211 Geneva 23, Switzerland 
8 Department of Physics, University of Manchester, Manchester, M60 IAD, United Kingdom 
9 Centre de Spectrometrie Nucleaire et de Spectrometry de Masse, F-91405 Orsay, France 
10 Institut de Physique Nucleaire d'Orsay, F-91406 Orsay, France 

11 GSI Helmholtzzentrum fur Schwerionenforschung, D-64291 Darmstadt, Germany 

12 Department of Physics, P.O. Box 3000, Fin-90014 University of Oulu, Finland 

13 Institut Laue-Langevm, F-38042 Grenoble, France 
14 Oliver Lodge Laboratory, University of Liverpool, Liverpool, L69 7ZE, United Kingdom 
15 KTH-Royal Institute of Technology, SE- 10044 Stockholm, Sweden 
16 Institute of Physics, Slovak Academy of Sciences, Bratislava 845 11, Slovakia 
17 Joint Institute of Nuclear Research, 141980 Dubna, Moscow Region, Russia 
18 Centre d'Etudes Nucleaires de Bordeaux Gradignan, F -33115 Gradignan, France and 
19 Service de Physique Nucleaire Theorique, Universite Libre de Bruxelles, B-1050 Bruxelles, Belgium 

(Dated: October 12, 2010) 

In-source resonant ionization laser spectroscopy of the even-^4 polonium isotopes 192 - 210 . 216 > 218 p o 
has been performed using the 6p 3 7s 5 S , 2 to 6p 3 7p 5 P2 (A = 843.38 nm) transition in the polo- 
nium atom (Po-I) at the CERN ISOLDE facility. The comparison of the measured isotope shifts 
in 20 °- 210 p o w itn a previous data set allows to test for the first time recent large-scale atomic cal- 
culations that are essential to extract the changes in the mean-square charge radius of the atomic 
nucleus. When going to lighter masses, a surprisingly large and early departure from sphericity is 
observed, which is only partly reproduced by Beyond Mean Field calculations. 

PACS numbers: 21.10.Ft, 27.80.+W, 29.38.-c, 31.15.-p 



In spite of substantial progress in understanding the 
structure of light atomic nuclei using ab-initio calcula- 
tions [l[ , heavy nuclei are still modeled using shell- model 
calculations, symmetry-based approaches or mean-field 
descriptions. These concepts address an interesting 
aspect of the finite ensemble of strongly interacting 
fermions that makes the atomic nucleus, namely the sub- 
tle interplay between microscopic (individual nucleon) 
and macroscopic properties (collective behavior). The 
lead isotopes are a good example: while the doubly 
closed-shell nucleus 208 Pb, with Z = 82 and N = 126, is 
a textbook case of a shell-model nucleus 0, signatures 
of particle- hole symmetries have been found [3j. Fur- 
thermore, shape coexistence has been observed in the 
neutron-deficient isotopes, whereby states with different 
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shapes occur at low excitation energy [4| . In spite of the 
latter effect, recent studies have shown that the ground 
state of the very neutron-deficient lead isotopes remains 
essentially spherical, confirming the robustness of the 
Z = 82 shell closure in those isotopes In the neigh- 

boring mercury isotopes (Z = 80), a significant deviation 
from sphericity has been deduced around N = 104, mid 
way between N — 82 and N = 126, from, e.g., the mea- 
surement of the changes in the mean-square charge radius 
(6{r 2 )) [7J. From symmetry arguments a similar effect is 
expected to occur in the neutron-deficient polonium iso- 
topes (Z = 84) d,©]. 

By approaching N — 104, the intrusion of a presumed 
oblate band was seen at low energy in 192 Po [HI EH- In 
190 Po, evidence for a prolate configuration in the polo- 
nium isotopes was obtained 12j , which is believed to be- 



come the ground state in 188 Po [ljj [14|. In order to un- 
derstand the change in the configuration, a direct mea- 
surement of ground-state properties independent of nu- 
clear models is crucial. The <5(r 2 } can be extracted from 
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isotope shifts measurements using laser spectroscopy in a 
model-independent procedure, provided that the atomic 
parameters (electronic factor F and specific mass shift 
Msms) are known [l^]. However, since there is no stable 
isotope of polonium, those measurements are extremely 
challenging. In particular, the information on the atomic 
structure of the polonium atom itself was hitherto rather 
scarce and limited in precision [lH ]. Hence a previous 



study of the S(r 2 ) in ' 



'Po by laser spectroscopy [Tj 



had to rely on a nuclear model (Finite Range Droplet 
Model - FRDM [l8|) in order to extract the 8(r 2 ) from 
the isotope shifts. 

In the last 20 years, the progress in atomic structure 
theory has enabled large-scale calculations for open-shell 
atoms and ions 19] providing the necessary input to an- 
alyze laser-spectroscopic data. However, one still has to 
assess the accuracy and reliability of those calculations, 
especially in heavy systems, as they are needed to access 
the atomic and chemical properties of the super heavy 
elements. 

In this Letter, we report on the measurement of the 
isotope shifts of the neutron-deficient, even- A polonium 
isotopes from 210 Po down to 192 Po (2~i/ 2 = 33 ms) and 
of the neutron-rich isotopes 216 > 218 Po using in-source res- 
onant photoionization spectroscopy. Combining the iso- 
tope shifts of the present and previous data set [13] al- 
lowed us to test large-scale atomic calculations for the 
first time and to evaluate their accuracy. The measured 
shifts and calculated atomic parameters are then used to 
extract S(r 2 } and those are discussed in terms of nuclear 
models. 

The polonium isotopes were produced at the CERN 
ISOLDE facility by 1.4 GeV proton-induced spallation 
reactions of a thick depleted \JC X target as described 
in Ref. [20j j over two experimental campaigns using the 
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194- 



204p o m 

For A > 204, 



resonant ionization laser ion source 

2007 and 19 2 ,196,206-210,216,218p o in 20 09 

the surface-ionized francium isotopes overwhelmed the 
beam. The A > 204 polonium isotopes were instead 
studied without proton irradiation, using precursors pre- 
viously accumulated in the target. Without proton irra- 
diation, the short-lived, fast-releasing francium contami- 
nation disappears as it does not have any precursor. The 



isotopes 



206-210 



Po were obtained via the /3 + /EC decay 



of the isobaric astatine nuclei while the isotopes 216 > 218 Po 
were produced via the a decay of 224 Ra and 222 Rn, re- 
spectively. 

The a-decaying isotopes 19 2-i96,2i6,2i8p were stud- 
ied as described in Ref. [2pj ] . The /3-decaying isotopes 
2O0-204p Q were gent tQ the ISOLDE tape station and 

their decay measured using a plastic scintillator and a 
single coaxial HPGe detector. The isotopes 206 - 210 p o 
were directly counted in a Faraday cup. The measure- 
ments were performed by scanning the narrowband laser 
at 843.38 nm from the 6p 3 7s 5 S 2 to the 6p 3 7p 5 P 2 ex- 
cited state of the ionization scheme plj and monitoring 
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FIG. 1: From top to bottom: example of a laser scan around 
the 843.38 nm transition between the 6p 3 7s 5 52 and 6p 3 7p 
5 P2 atomic excited states in even-A Po-I for A = 218 down 
to A = 192. For A = 192, several scans have been added 
together to obtain sufficient statistics. The line is the result 
from a fit of those spectra using an asymmetric Voigt profile; 
more information on this asymmetry can be found in the text. 



the yields as a function of the frequency. The frequency 
scans are shown in Fig. [TJ The shape of the resonance 
can be described as a deformed Voigt profile. The 2007 
data were analyzed as described in Ref. [f| . For the 2009 
data, an asymmetry in the fit function was introduced 
through a different Lorcntzian width parameter on each 
side of the resonance. The difference comes from the use 
of different pump lasers between the two runs. The de- 

are 



duced isotope shifts with respect to 210 Po (Sv^: 210 ) 
presented in Table fl] 

The large overlap in the studied isotopes between the 
present data set and the previous study using the 255.8 
nm transition FjJ allows the formalism of King to be used 



to study the atomic properties of polonium by comparing 
modified isotope shifts (see Fig. [5]) [lfij]. The slope of this 
so-called King plot is the ratio of the two F while the y 
intercept is a linear combination of the Msms from each 
transition. The Msms and F can also be calculated on 
the basis of the Dirac-Coulomb-Breit Hamiltonian and a 
Fermi-like distribution of the different isotopes. A series 
of relativistic configuration interaction calculations have 
been carried out with systematically enlarged wave func- 
tion expansions within a restricted active space, includ- 
ing the polarization of the electronic core and single and 
double excitations into three additional layers of correla- 
tion orbitals (n = 8, 9, 10). Reasonable convergence with 
the size of the wave functions was obtained, especially 
for F, while the Msms values appear more sensitive to 
the details of the calculations. The results are shown in 



TABLE I: Isotope shifts 6u^p 10 in the 843.38 nm transition of 
Po-I and changes in the mean-square charge radii S(r 2 )f^ 210 
of the polonium isotopes with respect to 210 Po from this work. 
The systematic uncertainty on the 5{r 2 )^p 10 originates from 
the Specific Mass Shift. 
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FIG. 2: (Color online) King plot between the transitions at 
255.8 nm ([13, x axis) and at 843.83 nm (present work, y 
axis) for 20 °- 210 p o . The solid line is a linear fit through the 
data points; the dotted red line is the calculated relation from 
the large-scale atomic calculation and lies within 1 a of the 
fit. 



Table HB The red dotted line in Fig. [5] displays the corre- 
sponding relation according to the calculated parameters 
and lies within 1 a of the fitted trend. The good agree- 
ment between the slope of the calculated parameters and 
that of the fit shows the predictive power of the calcula- 
tions for F. The difference in the y intercept, however, 
raises some questions on the theoretical accuracy of the 
calculated Msms- 

The 5{r 2 ) of the even-A isotopes 192 - 210 > 2i6,2i8p Q were 
extracted using those parameters and a 0.932 correction 
for higher moments [22J. In order to take into account 
the uncertainty of the different Msms m the calculation, 
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TABLE II: Calculated atomic electronic factors F and specific 
mass shifts Msms- 
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FIG. 3: (Color online) (top) 5(r 2 } for the even-i? isotopes 
from platinum (Z = 78) to radium (Z = 88) [1, 0, H34I11 ■ The 
solid black line represents the predictions from the spherical 
FRDM [2 using the second parametrization from Ref. [23T ]. 
(bottom) Difference between the measured 5{r 2 } to the spheri- 
cal FRDM. The dotted lines represent the Beyond Mean Field 
calculations with the SLy4 and SLy4* (with reduced pairing) 
interactions pffll. [30l|. 



a systematic error was introduced. It was deduced as the 
difference between the 5(r 2 ) values using only the calcu- 
lated atomic parameters for the 843.38 nm line and those 
obtained via the King plot and the calculated atomic pa- 
rameters from the 255.8 nm transition. The 8(r 2 ) (sec 
Table [l| are compared with the predictions from the 
spherical FRDM [l8j using the second parametrization 
from Ref. [23| (see Fig. [3]). On the neutron-deficient side, 
a surprisingly large deviation from sphericity can be seen 
starting from liijPo that becomes increasingly marked 
for the lighter isotopes. The deviation is larger in magni- 
tude and occurs for larger neutron numbers than in the 
Z < 82 isotones. The data in the neutron-deficient radon 
and radium isotopes [24T ] do not extend far enough in the 
neutron-deficient side to compare with the polonium iso- 
topes. 

In order to understand the unexpectedly large and 
early deviation from sphericity in the polonium isotopes, 
the S(r 2 ) have been calculated using the same Beyond 
Mean Field method as in Ref. [i^, \3(^. The most im- 
portant feature of the method for this study is that the 
ground-state wave function is constructed as a superposi- 
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tion of mean-field wave functions corresponding to a large 
set of axial quadrupole deformations, projected on angu- 
lar momentum and particle number. The coefficients of 
the expansion are determined by varying the energy cor- 
responding to a Skyrme energy density functional. The 
SLy4 Skyrme parametrization has been tested together 
with the effect of a reduced pairing strength (SLy4*). 
Within this framework, one cannot assign an intrinsic 
deformation to the wave functions. Instead, they are a 
mixture of mean-field states of different deformation and 
therefore different radii. In general, deformed configu- 
rations have larger radii than spherical ones. The two 
main effects that increase the radii of neutron-deficient 
polonium isotopes, compared to the global A 1 ^ trend 
set by spherical configurations, are the spread of the 
collective wave function in deformation space and the 
shift of the dominant configurations from near-spherical 
to oblate. The increasing softness of the deformation 
energy surfaces, when going down from 210 Po to 194 Po, 
leads to collective ground-state wave functions of increas- 
ing spread, but which remain centered around spherical 
shapes. For 192 > 190 p Oj the ground-state wave function 
becomes centered around an oblate minimum in the de- 
formation energy surface and the contribution from near- 
spherical configurations (or smaller radii) becomes sup- 
pressed. 

The calculated 5 (r 2 } are compared with the experimen- 
tal data at the bottom of Fig. [31 after subtraction of the 
FRDM value. There is a qualitative agreement between 
theory and experiment, although deficiencies exist, es- 
pecially for 192 . 194 p 0; where the data indicate a stronger 
deviation from sphericity. The effect of a reduced pairing 
strength, which clearly improves the agreement between 
theory and experiment for the lightest nuclei, is meant 
to put a larger weight on deformed oblate configurations. 
One still needs to construct more flexible energy func- 
tionals to correct the deficiencies of the actual ones. 

Finally, the neutron-rich isotopes 216 < 218 p show a 
clear break from the trend of the polonium isotopes below 
N = 126. The magnitude of this kink is similar to what is 
observed in the neutron-rich neighboring lead [Z = 82) 
[jjj ]. bismuth (Z = 83) 31] and heavier isotopes. The 
underlying mechanism is still an open question. 

In conclusion, in-source resonant ionization laser spec- 
troscopy has been performed on the polonium isotopes 
from 192 Po to 218 Po. The overlap with the previous data 
set available in the literature has allowed testing of the 
large-scale atomic calculations for the F and the Msms, 
although the latter remain slightly less well determined 
for such heavy atoms. This first experimental evidence of 
the reliability of such calculations is crucial for the study 
of laser spectroscopy in complex systems [la ]. Moreover, 
it shows that reliable information can be extracted for the 
very heavy elements, where limited or no atomic data are 
available yet. 

The 5(r 2 } of the even- A polonium isotopes 



192-210, 216, 2i8p nave Dee n extracted and compared 
with systematics of this region and recent calculations. 
An unexpectedly large departure from sphericity was 
observed compared with the Z < 82 isotones. Com- 
parison to Beyond Mean Field calculations indicate 
that the coexistence of the different shapes at low 
excitation energies leads to a very soft nature of the 
most neutron-deficient polonium nuclei. The different 
trend with respect to the Z < 82 nuclei might suggest 
that high-j orbitals occupied by the protons play a 
critical role. The study of the more neutron-deficient 
radon [Z = 86) and radium [Z = 88) isotopes could 
eventually shed more light on the aspect. 
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